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The Laguerre-Gauss modes are a class of fundamental and well-studied optical fields. These stable, 
shape-invariant photons - exhibiting circular-cylindrical symmetry - are familiar from laser optics, 
micro-mechanical manipulation, quantum optics, communication, and foundational studies in both 
classical optics and quantum physics. They are characterized, chiefly, by two modes numbers: the 
azimuthal index indicating the orbital angular momentum of the beam - which itself has spawned 
a burgeoning and vibrant sub-held - and the radial index, which up until recently, has largely been 
ignored. In this manuscript we develop a differential operator formalism for dealing with the radial 
modes in both the position and momentum representations, and - more importantly - give for the 
hrst time the meaning of this quantum number in terms of a well-dehned physical parameter: the 
“intrinsic hyperbolic momentum charge”. 


I. INTRODUCTION 

The ability to shape optical beams has increased dra¬ 
matically in recent years thanks to a number of newly 
developed tools and techniques. This progress has been 
driven by the many promising potential applications of 
the orbital angular momentum degree of freedom, as well 
as more generalized beam shaping. In 1992 Allen et 
al. [T] showed the connection between the azimuthal in¬ 
dex of a Laguerre-Gauss (LG) mode (the natural beam¬ 
like light modes exhibiting circular-cylindrical symme¬ 
try) and the physical quantity of orbital angular momen¬ 
tum (0AM). This physical parameter is useful for doing 
micro-mechanical work, is a well-defined quantum num¬ 
ber which applications to foundations, and informatics, 
and also is important in the foundational study of optics 
(see Ref. [5] and references therein). 

However, despite the vast amount of attention 0AM 
has received over the past two decades very little re¬ 
search has been conducted on the other mode number 
of LG photons - the radial index. Much of the time it 
is regarded as little more than unwanted noise which un¬ 
fortunately arises when one is trying to produce “pure” 
0AM modes. In a previous note we referred to this as 
the “forgotten quantum number” [3]. In this paper we 
will summarize the previous work by us and by others 
HIS], then significantly extend and simplify these pre¬ 
vious analyses, elucidate the physical “meaning” of this 
quantum number, and give a brief prospectus for the use 
of the radial number in future quantum technologies. 

Laguerre-Gauss beams are typically found by tak¬ 
ing the paraxial wave equation and finding solutions in 
the circular-cylindrical coordinate system. The paraxial 
equation is the result of making the assumption that the 


‘Electronic address: billplick@gmail.com 


beam is not highly divergent or focused. Ironically, we 
find the mathematical analysis of the radial number sim¬ 
plifies signihcantly when this assumption is removed and 
exact solutions of Maxwell’s equations are considered. 
This is done by taking the photonic wave-function in 
momentum coordinates as derived through the Riemann- 
Silberstein vector formalism mi- However, on the other 
hand, the physical meaning does remain clearer in the 
position representation, an oddity we will discuss. 


Our conclusion is that the radial index of Laguerre- 
Gauss photons is a compound physical parameter influ¬ 
enced by - but not influencing - other fundamental pa¬ 
rameters of the mode, as well as an additional property 
intrinsic to itself: the hyperbolic momentum; which is 
a kind of mathematically well-formed “radial-like” mo¬ 
mentum with subtle and interesting properties. Thus we 
call the radial index the “hyperbolic momentum charge”. 
The hyperbolic momentum is the result of the restric¬ 
tion that the radial coordinate is only defined for values 
greater than zero, as we will explain. 


This paper is organized as follows: In the following 
section we review the LG beams, in section 3 we derive 
the radial mode operator in the paraxial, position-space 
representation and introduce the hyperbolic momentum 
operator. In section 4 we briefly summarize the back¬ 
ground for the exact quantum momentum-space wave- 
function of the photon, namely, the Riemann-Silberstein 
vector - a more detailed derivation being included in the 
appendix. In section 5 we derive the momentum-space 
formalism for the radial modes. In section 6 is the raison 
d’etre of this manuscript: the physical interpretation of 
the preceding mathematics as well as a discussion of some 
potential applications of the radial quantum number. In 
section 7 we summarize and conclude. 
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II. LAGUERRE-GAUSS BEAMS AND THE 
ORBITAL ANGULAR MOMENTUM OF LIGHT 

The equation of a LG beam, under the paraxial as¬ 
sumption, in circular-cylindrical coordinates (r, (j), z) is, 
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Where n and I are the radial and orbital angular mo¬ 
mentum quantum numbers, respectively; L^n is the gen¬ 
eralized Laguerre polynomial of order n and degree \l\. 
The functions R^, and ipg are the beam waist, radius 
of curvature, and Gouy phase of the fundamental beam, 
and are given by 
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Where k = 2'K jX is the overall wavenumber of the beam, 
and Wo is the beam waist at z = 0 (defined where the 
beam is narrowest). The beam or photon is then com¬ 
pletely defined by four numbers n, I, k, Wo- In Figure 
[2 we plot some LG beams for various parameters. The 
LG beams are solutions to the paraxial wave equation 
(PWE) 


The 0AM is also sometimes called the “winding num¬ 
ber” as it increments the number of helical equal-phase 
surfaces of the photon (see FigQ. If a particle absorbs a 
photon with both spin (arising from the polarization of 
light) and orbital angular momentum the spin will cause 
the particle to rotate about it’s own axis, whereas the 
0AM will cause it to rotate about the optical axis of the 
beam m- 


HI. THE RADIAL MODES - OPERATOR 
FORMALISM IN THE PARAXIAL REGIME 

The orbital angular momentum of Laguerre-Gauss 
beams has received an extensive amount of attention, 
however the radial index [n in our notation) has been 
the subject of only a handful of manuscripts. The ra¬ 
dial index is so-called because the intensity pattern of 
LG beams exhibit n -I- 1 concentric rings if Z 0 (there 
are n rings around a central dot for 1 = 0). Also, the 
phase structure of the beam displays n concentric radial 
discontinuities with no smooth transitions at z = 0 (see 
FigQ. This is in contrast to the OAM/azimuthal coor¬ 
dinate where the phase does have smooth transitions at 
z = 0. Unlike the 0AM operator, Eqj^ the differential 
operator for the radial number has not been studied until 
very recently. 

The operator formalism for the radial modes of LG 
beams was first derived by Karimi and Santamato [4] in 
2012 via group-theoretic techniques, and then indepen¬ 
dently by Plick et al. in 2013 [3] via the method described 
below. 

We start with the Laguerre polynomial. There ex¬ 
ist a series of relations between Laguerre polynomials 
of varying order and degree (the so-called “three-point 
relations”). One such relation is 


d 

VlE -2ik—E = 0. (5) 
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Where Vj is the transverse Laplacian and A is a com¬ 
plex electric scalar field (i.e. we assume the electric field 
vector points in the same direction at every point in the 
transverse plane). The paraxial wave equation is a ver¬ 
sion of the full wave-like Maxwell’s equations under the 
small angle approximation for propagation. In this paper 
we will consider both the cases where this approximations 
is met, and when it may not be. 

In the paraxial limit the angular momentum of light 
separates out into a spin and orbital part, both of which 
are well-defined m- The form of the differential 0AM 
operator in the paraxial limit, about the direction of 
beam-propagation, is well-known and straightforward 
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Which, when combined with the rule for differentiation 
of the polynomials. 
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Given this differential relation, it is possible to arrive 
at a relationship between the full LG modes by left- 
multiplying the other factors (the non-polynomial terms) 
in the LG function Eq. 0 onto Eq. 0, and commuting 
those factors past the differentials on the left. Doing this 
we obtain 
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FIG. 1: The transverse spatial profiles of nine different Laguerre-Gauss beams in both intensity (left) and phase (right) at 
z = 0, for three different values of the radial index and three different values of the azimuthal index. The diameter of all the 
plots is 6mm. 
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Now if we make the the identification x = 2r^ and 
perform a change of coordinates, we obtain 
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Where, iVoLG„/(r, (/), 0) = nLG„i(r, 0), and we have 
defined the differential n-mode operator for z = 0. A 
trivial transformation shows that this is indeed the oper¬ 
ator derived in [4]. 

Further simplification is possible by identification with 
other, better known, operators yielding the opportunity 
for some physical insight 





( 12 ) 


The operator is composed of the transverse Laplacian in 
cylindrical coordinates, and thus has a direct relation¬ 
ship to the transverse momentum squared (up to some 
constants). For completeness this Laplacian is defined as 
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Equation (12) also contains the 0AM operator (with 
h = 1), as well as a term reminiscent of a harmonic po¬ 
tential. It commutes with the 0AM operator, but not 
with the Hamiltonian for free propagation due to the fi¬ 
nal potential-like term. This indicates that n is not a 
conserved quantity of the beam with respect to evolution 
through free space. This point will be discussed further 
in section 6. 

It is also useful to derive a general operator for any 
value of z. Gontinuing with the same methods as above 
we have 


iz d Lz 1 / r“ , , 


The differences being that the pre-factor of Laplacian 
is now the z-dependent beam waist and an additional 
term resultant from the phase imparted from the radius 
of curvature, (this can be seen from the fact that if the 
LG beam is written without the radius of curvature phase 
factor then this term does not appear). Interestingly, the 
final term is unchanged and remains z-independent. 

The z-dependent form in Eq.(14) offers another insight. 


The second term (radius-of-curvature term) can be iden¬ 
tified as the operator for a well-formed radial momentum 
known as the “hyperbolic momentum”. 















4 


It has been known for a long time [H] that the operator 
for a strictly radial momentum can not be well defined. 
To see this consider the most direct attempt: = —iHdr, 

and its action on the radial coordinate operator: 


( 15 ) 

But the domain of f is only the non-negative real num¬ 
bers. The operator can take r out of this domain 
and is thus not well-formed. An operator which can be 
well-defined in the circular-cylindrical coordinate system 
is the hyperbolic momentum 


P}j = — ih. 



(16) 


Which, up to constants (and the z coordinate) is the 
second term in Eq.(14). The action of this on the radial 
coordinate is 
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Thus the hyperbolic momentum generates dilations, not 
linear translations, and can not cause the radial coor¬ 
dinate to be negative no matter the value of 7 . As 
linear momentum is associated with invariance under 
translation, hyperbolic momentum is associated with in¬ 
variance under scale transformations. However, since 




—2ihS/1 hyperbolic momentum is not a con¬ 


served property of paraxial photon propagation (unlike, 
obviously, the linear-z momentum and the 0AM). It is 
clear that in our case the position r can not be negative, 
that is r G . A convenient method of dealing with this 
is to make the transformation rj = ln(r), which is only 
dehned for positive positions. The conjugate variable to 
rj is the hyperbolic momentum. With this construction 
it is unsurprising that the hyperbolic momentum should 
appear in our formalism. The hyperbolic momentum has 
many other interesting properties, for a detailed inves¬ 
tigation the interested reader is referred to Ref. m and 
Ref.[T4]. 

For the case of a LG-beam the expectation value of the 
hyperbolic momentum increases linearly as a function of 
the propagation distance z. Thus the expectation value 
of the second term in Eq|^ as a whole has a quadratic 
scaling in z. The hyperbolic momentum is always zero 
where the beam-waist is narrowest, since at this point the 
beam is neither dilating nor contracting (see Figj^. If 
we look at the same quantity as a function of beam waist 
Wo we see an exponential decay. This later effect can be 
understood from the fact that as beam waist increases so 
does the degree of collimation of the beam, thus the beam 
dilates less (see Figj^. These functions can be directly 
calculated as 


(Ph) 
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FIG. 2: The expectation value of the hyperbolic momentum 
as a function of propagation distance for five different val¬ 
ues of the radial index. In all cases the scaling is linear, 
passing through the origin where the beam-waist is small¬ 
est. Note that this does not include the pre-factor of —zfkWo 
in Eq.( |14[ |. For that whole term together the z-scaling is 
quadratic. Which, in the far field is equivalent to the scaling 
of the Laplacian term. 



FIG. 3: The expectation value of the hyperbolic momentum 
as a function of beam-waist at origin for five different values of 
the radial index. The hyperbolic momentum has a finite value 
at Wo (not pictured due to scaling in order to show contrast 
between modes). 


See Figures and 1^ for an illustration of the numerical 
results calculated from Eq.(18|. 

We will examine in more detail the physical interpre¬ 
tation of the operator in Eq.(14) in section 6 . 

It has been previously noted by Karimi and Santam- 
ato in Ref. [4] that the 0AM and the radial index are 
inextricably linked. From the referenced paper: “0AM 
and radial intensity distribution are strictly correlated, 
and different 0AM generators produce specific (and dif- 
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ferent) distributions of radial modes.” (We will comment 
again on this fact in section 6). In that paper “radial co¬ 
herent” (displaced vacuum) and “intelligent” (minimum 
uncertainty) beams are also derived. The latter of which 
can be generalized to squeezed states which is done by 
Karimi et al. in Ref. [5] via a sophisticated algebraic tech¬ 
nique employing the raising and lowering operators on 
the radial number - which form a lie algebra. It is also 
possible to alternatively define the coherent state as the 
eigenstate of the lowering operator, which is also done in 
Ref. [3]. Unlike light in the number basis - where all three 
definitions are connected via the familiar coherent state 
I a) and its generalization the single-mode squeezed state 
- in the radial representation, all three concepts result in 
distinct beams. 

There has also been experimental work done on the 
value of the radial index as a quantum number. Again 
this was carried out by Karimi and collaborators in 
Ref. [5]. They showed via Hong-Ou-Mandel interference 
that the radial degree of freedom is indeed a quantum 
number and could in principle be used in quantum- 
information tasks. However, a small caveat is that care 
must be taken to ensure that the beam waists are equiv¬ 
alent. From that paper: “The chosen basis is beam-waist 
dependent; an eigenstate for a specific beam waist turns 
into a superposition of radial modes for any other beam 
waist.” 

Several experiments have taken advantage of both the 
0AM and radial indices, for example Ref. [15]. It has 
also been shown |16j that via the use of both quantum 
numbers, entangled states of very high dimension may 
be produced. Even as high as 103 dimensions. Research 
has also been done on efficiently producing radial modes 
[HI- It has also recently been shown that LG beams of 
high radial indices exhibit self-healing properties in the 
same way that Bessel-beams do [18]. 


IV. THE RIEMANN-SILBERSTEIN VECTOR, 
EXACT SOLUTIONS TO MAXWELL’S 
EQUATIONS, AND THE PHOTONIC 
WAVE-FUNCTION 


to Ref. [7] and Ref. [5]. In the appendix we review this 
derivation in the interest of completeness and notational 
consistency. The previously acquainted reader may skip 
it. 

It is convenient to work with the scalar function x 
(which contains the full information of F, up to choice of 
gauge, as explained in the appendix) as opposed to the 
full RS vector F. Given the formalism from the appendix 
we can directly write a general beam in the Bessel-Gauss 
basis as 


xl,{rA,z,t)= J ( 20 ) 


Using Eqs.( 51|M ). Note that we take a different defini¬ 
tion than is used in Ref. [S] as we pull the azimuthal phase 
and Jacobian into tp. 

Now, derivation of the Laguerre-Gauss modes is fa¬ 
cilitated by a change of momentum coordinates k± = 
{k ± kz)/2 and also time t± = t ± zjc, yielding 


Xm(U </’,z,t) 



xJ™(2rVA:+fc_). (21) 


It is of key importance to note that since y is ef¬ 
fectively the position-space wave-function (given the 
caveats described in the appendix) that ip'^{k+,k-) is 
the momentum-space wave-function in the beam-like (i.e. 
Bessel) basis. Now, if we make the choice 


r{.k+,k.,k^) = 


j n+^ _ 

xk ^ e 


■k. 


( 22 ) 


Where n is some integer, and w is some real number, we 
can solve Eq.(21) exactly: 


Here we briefly outline some previously existing math¬ 
ematical formalism we will need. A more complete ver¬ 
sion is included in an appendix. The Riemann-Silberstein 
(RS) vector is a representation of the electromagnetic 
field given by 


Xm(U (p,Z,t) 
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F = ^ [e + icB'j . (19) 

Where c is the speed of light in vacuum, eo is the per¬ 
mittivity of vacuum and E and B are the electric and 
magnetic fields, respectively. The vector consists of three 
complex numbers and many useful quantities can be cal¬ 
culated directly with simple equations. For a full re¬ 
view on the RS vectors the interested reader is directed 


Which is clearly a form of the Laguerre-Gauss beams, 
with Af a normalization constant, m an angular mo¬ 
mentum number (we avoid calling this the orbital an¬ 
gular momentum since the 0AM is only completely well- 
defined in the paraxial limit), n the radial number, and 
a(t+) = -\- iac^t+IQ, - a form of the complex beam 

parameter. Thus we can identify Eq.([2^ as the ex¬ 
act momentum-space wave-function of a Laguerre-Gauss 
beam (up to normalization constants). 
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V. THE MOMENTUM-SPACE FORMALISM 
FOR THE RADIAL MODES 


In the previous section, and in the appendix, we sum¬ 
marized the mathematical tools from Ref.[S] which we 
will need for the investigation of the radial number as it 
exists for the exact photon wave-function. 

From the momentum representation of Laguerre-Gauss 
beams in Eq. (22 1 , we can easily derive an radial momen¬ 
tum operator. 


As well as k = « k^- Now, again take Eq.(22) 

- or use directly Eq. (291 - and re-write it in terms of 
these more familiar momentum coordinates 


r 








K (32) 


And making use of the delta function, we have up to 
constants. 
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Therefore, one finds the (general, non-paraxial) radial 
momentum operator as 



i d 
2a dk^ 



(26) 


In order to interpret the operator, one can transform it 
to polar coordinates. By using the relations 


r{kuk,,k^) oc (33) 

From this, unlike in the position-space paraxial-regime 
case, it is quite straightforward to find the operator which 
returns the value n. By simple inspection we can write 


N'k 



i d 
a dk^ 


-H w'^k; 


(34) 


Which has the most straightforward form. The first term 
is similar to the hyperbolic momentum operator, however 
it lacks the i in front. This is non-trivial as without the 
i it’s action on the radial coordinate is given by 


k - (27) 

t* = (28) 

We find directly the polar-coordinate representation of 
the radial-momentum operator in momentum space: 



(fc - kz){ 


d 

dkz 


1 

k 



(29) 


This is the most exact form of the operator. Now, if we 
take the case of paraxial photons (which is in most cases 
the situation of interest), the momentum along the di¬ 
rection of propagation will be many orders of magnitude 
larger than the transverse momentum. So if we Taylor- 
series expand the square-root in k- and k^ about k^ we 
have 



fcz. (31) 


Where the radial position coordinate is defined in the 
momentum representation as 


fs-ttA. (36) 

Obviously this is problematic as the radial coordinate 
must remain a real number. Also, without the i the the 
first term (by itself) is not in general hermitian. It’s 
potentially interesting to note that if the action ( 7 ) is 
a multiple of tt then it generates a real (periodic) scal¬ 
ing. Whether this can have some physical meaning is 
unknown to us as of yet. 

The second term is again the 0AM operator, scaled by 
the helicity. The final term is the transverse laplacian in 
momentum coordinates scaled by the beam waist at ori¬ 
gin. There are also a couple other advantages of writing 
the operator in this way. It is propagation invariant, that 
is, this is the operator for all time (as the momenta are 
conserved quantities). Also, the transverse momentum 
coordinate may be rescaled (as k[ = wkt) io make the 
operator independent of the beam waist - should this be 
desirable. 

We want to test is its hermiticity (in total). Does 
(iV^' 0 '^|'!/'‘’’) = ('!/)‘^|7V(.^/:‘’’)? a straightforward computa¬ 
tion shows that if we desire real eigenvalues, then N'j^ 
can not, in general, be hermitian. It is hermitian in¬ 
stead only on a restricted class of wavefunctions where 
tp{kt,k^) = Q{kt)^{k^), where Q{kt) is a real function. 
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Only in such beams can Ni~ be considered a valid observ¬ 
able. Some examples of such beams are the Bessel beams 
and the Laguerre-Gauss beams. Which is our main in¬ 
terest. 


VI. PHYSICAL INTERPRETATION AND 
TECHNOLOGICAL PROSPECTS 


Now that we have laid down the mathematical foun¬ 
dations for the radial index operator we can inves¬ 
tigate “meaning” of the radial index. Though the 
RSV momentum-representation version of the operator, 
Eq.(34|, is accurate and simple, it’s interpretation is 


more difficult to parse than the paraxial coordinate- 


representation operator Eq. (14). This is indicative of the 


fact that the radial-index operator is not representative 
of a general quantum observable like it’s “partner” the 
0AM operator. Instead n can be considered a proper 
observable in a restricted class of scenarios. More specif¬ 
ically for beams whose momentum-space wave-function is 
separable into radial and angular portions, and for which 
the radial part is always real. The radial-index opera¬ 
tor also requires additional information to write down 
- specifically the beam-waist at origin, the propagation 
distance, and the orbital angular momentum. 

This latter point has implications for future potential 
uses of the radial index as a tool in quantum technologies. 
Suppose we wish to measure the (unknown) radial index 
of a single photon. A common technique for doing this 
would be to shine the photon on a phase pattern which 
takes a photon of the radial index in question and turns it 
into the fundamental Gaussian mode, this photon then 
either couples - or does not couple - to a single mode 
fiber. This constitutes a projective measurement onto 
the chosen radial number. However, in order to choose 
the correct phase pattern, the 0AM index, the beam 
waist, and the distance of propagation must be known. 
Note that this is not merely a problem with a particu¬ 
lar experimental implementation of radial-index-sorting, 
but a general problem. If we examine the paraxial version 
of the operator - Eq. (14) - we can see that it contains 


all three of the other parameters {z, Wq, and 1) mean¬ 
ing that we must know these values a priori (or measure 
them simultaneously) in order to get a meaningful result 
from a radial-index measurement. It may also be pos¬ 
sible to perform a measurement if the other parameters 
can be measured passively, that is, non-destructively for 
the same photon (for example, 0AM may be passively 
measured with a mode sorter M)- 

To further illuminate this problem it should be noted 
that Laguerre-Gauss beams of different radial-index do 
not have zero overlap if they are at different distances in 
their propagation. Mathematically, that is 


J rdr J d(j)LGli{r,(j),z)LGn'i{r,(j),z') ^0. (37) 



Propagation distance mismatch : 2 (m) 


FIG. 4: Graphs of the overlap function of two Laguerre-Gauss 
modes with un-equal radial indices as a function of the mis¬ 
match in the distance the modes have propagated. In both 
graphs the black line represents the mode’s overlap with it¬ 
self. In (a) we show the n = 0 mode’s overlap with the five 
next higher indexed modes. It is clear from the graph that 
very quickly the depicted modes no longer add to one. This 
is due to the fact that as the beam propagates it acquires 
non-trivial overlaps with more and more modes. After only 
about 15 meters 20 modes are necessary to describe the beam 
to good accuracy (not shown). In (b) we show the overlap 
of an n = 5 mode with a mode of higher index and a mode 
of lower index. Scaling is different than in (a) to highlight 
detail. The oscillatory behavior is due to the fact that after 
“moving” from one mode to an adjacent one, the photon may 
then “move back” as it continues to spread in the radial-index 
space. 


In Figurej^we plot the overlap of a beam of a particular 
n with several others for increasing propagation-distance 
mismatches. 

However if these problems can be overcome (which can 
be achieved by careful calibration) then indeed the radial- 
index can be used as a carrier of quantum information in 
realistic scenarios. This has been demonstrated recently 
by two separate experiments, as discussed earlier. 

These experiments show that the radial index can be a 
valuable resource in quantum experiments, but one needs 








to be careful to satisfy the subtle requirements in the 
measurements. 

The maxim “One man’s noise is another man’s data” 
may also hnd some application for the radial modes. To 
wildly speculate, it could be the case that the index’s 
sensitivity to distance and beam-width mismatch could 
find some application, and that by (for example) measur¬ 
ing one n-mode’s projection onto nil some information 
about propagation distance or dilation could be obtained. 

Returning to the question of physical meaning, it is 
more illuminating to consider the paraxial version of the 
operator, Eq.(14). Which we re-write below for the the 


convenience of the reader 


4 = + 


kw^ 


wf, 


(38) 


Where we have multiplied through by h and substi¬ 
tuted the second term with the hyperbolic momentum 
operator - Eq.(16|. There are four terms, three of which 


are related to fundamental properties of the beam which 
can be identihed with quantities other than - and inde¬ 
pendent of - the radial index. That is to say that these 
three parameters influence - but are not influenced by - 
the radial index. This leaves a single term which repre¬ 
sents the degree of freedom that the radial index truly 
represents. The first term in Eq.(38l is the transverse 


Laplacian scaled by the z-dependent beam waist and a 
numerical constant. In the free-space paraxial regime 
the transverse Laplacian alone is the Hamiltonian (and, 
of course, a constant of motion). The third term is the 
0AM. The hnal term produces - up to constants - the 
second moment of the radial position at z = 0, which is 
clearly related to the transverse spatial variance, i.e. the 
spatial confinement at origin. These are three indepen¬ 
dent parameters of the beam (and of individual photons 
which occupy these modes). 

It is the second term, however, which is most impor¬ 
tant for understanding the physical meaning of the ra¬ 
dial index. This is the hyperbolic momentum operator 
which generates dilations. As discussed previously it is a 
mathematically well-formed kind of “radial momentum 
operator”. It is however not a constant of motion (as 
is demonstrated in Eigj^, and thus does not have one 
of the quantities we would traditionally expect from a 
momenta. Nonetheless, if the propagation distance is ac¬ 
counted for in the operator itself (as is the case in iV^), 
then it is sensible to speak of it as an intrinsic prop¬ 
erty; here we see why the operator must be z-dependent 
in the coordinate representation and becomes somewhat 
unmanageable in the momentum representation. In other 
words, the radial index is best understood in the coor¬ 
dinate picture, as in that representation the fact that 
it is not conserved under time-evolution is easy to ac¬ 
count for, however the momentum representation deals 
best with quantities that are preserved for all time. 

The radially transverse momentum of a photon must 
be dependent on the total transverse momentum, as well 


as the orbital angular momentum (as the orbital part 
of the total transverse momentum needs to be compen¬ 
sated for); and on the degree of spatial confinement of the 
beam (because greater confinement must lead to higher 
divergence and greater radial momentum). However this 
does not account for intrinsic freedom of the hyperbolic 
momentum. The radial index can then be seen as repre¬ 
senting the promotion of the radial-like, hyperbolic mo¬ 
mentum itself to a quantum observable. We therefore 
believe that the radial index could best be described as 
the “intrinsic hyperbolic momentum charge”. 

We might have initially expected the radial index to be 
related to the radial momentum, as the azimuthal index 
is related to a rotational momentum, completing the set 
of two momenta necessary to move about in the trans¬ 
verse 2-D plane. We have found that, while Nz does 
not correspond to a full radial momentum operator (as 
no such well-formed operator exists), the hyperbolic mo¬ 
mentum operator - which is up to scalings and other ob¬ 
servables the operator we hnd - has a direct connection 
to the momentum along the radial coordinate. 

It has been shown m that the radial index is con¬ 
served in the process of spontaneous parametric down- 
conversion. More specihcally, in the ideal case, the two 
daughter photons will have the same radial index as each 
other. In the non-ideal case - hnite crystal length and 
width - there is a spread in the radial index of the down- 
converted photons. As we have seen in the research in 
this manuscript, the radial index can only be said to be 
well-dehned given a precise propagation distance, 0AM, 
and beam-width; therefore we gain insight as to why the 
effect of crystal length (and thus uncertainty in the point 
of origin of the down-conversion; and walk-off) should 
“blur” the resultant radial indices. This is in addition to 
a similar effect from a limited transverse spatial extent 
of the crystal - investigated in detail in Ref. [20] 

Since this paper and the previous work of ourselves and 
others have shown that both the orbital angular momen¬ 
tum and the intrinsic hyperbolic momentum charge are 
well-formed quantum numbers with a physical meaning 
(in the proper circumstances), it follows that that the 
mode indices of beams with different symmetries (non¬ 
circular cylindrical, e.g. elliptic or Cartesian) are also 
well-formed quantum numbers. These other mode indices 
are just algebraic combinations of I and n. Therefore, it 
should follow that these other numbers should have inter¬ 
esting physical meanings as well. Eor the case of parax¬ 
ial photons with elliptic symmetry - the so-called Ince- 
Gauss modes [25] - it has been shown that their mode 
indices can be be entangled [22] • Also, that their 0AM 
properties are non-trivial and show interesting complex¬ 
ity |2S]- It may be the case that these effects could be 
better understood by combining the concepts of orbital 
and hyperbolic momentum. 

There are still some unanswered questions - most sig¬ 
nificantly - it is not clear should be considered the con¬ 
jugate variable. Also, we do not yet understand why 
the “intrinsic hyperbolic momentum charge” (radial in- 
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dex) should take discreet values. Furthermore it is un¬ 
known why the hyperbolic momentum should be tied to 
the number of rings in the transverse spatial pattern of 
the beam. 

We would also like to briefly sketch some potential 
technological applications and concerns tied to the ra¬ 
dial index. As we have mentioned already it has been 
shown experimentally that the radial index is indeed a 
valid quantum number [^, and could be used for quan¬ 
tum communication. However the fact that different n- 
modes can have non-zero overlap if there is mismatch in 
the beam waist Wo, or propagation distance z, adds some 
additional complication, especially in realistic scenarios 
where complete control or knowledge of the distances and 
amount of focusing involved may not be available. How¬ 
ever, if we take the z = 0 form of the radial-index opera¬ 
tor, Eq.( 121, we find that this is exactly the Hamiltonian 
of a graded-index fiber (GRIN) [24]. GRINs are commer¬ 
cially available fibers that have a linear radial variation of 
their refractive indices. It is possible that a hyperbolic- 
momentum carrying photon could be matched with stan¬ 
dard optics to such a fiber (i.e. beam waist matched to 
the correct steepness of index variation) then the radial 
index would be intrinsically preserved under propagation 
in the fiber. This would allow for quantum or classical 
communication multiplexing using the radial degree of 
freedom in a widely-available kind of fiber. This is an 
exciting prospect that we believe merits further investi¬ 
gation. 


dial index but are influential on it. Put another way, for 
each value of these parameters (beam waist, orbital angu¬ 
lar momentum, propagation distance) there is a different 
representation of the radial-index operator. There is one 
remaining part of the operator which is not tied to other 
properties of the photon which we have identified: the hy¬ 
perbolic momentum which generates dilations. It is this 
term which corresponds to the true degree of freedom the 
radial index represents. We thus call the the radial index 
the “intrinsic hyperbolic momentum charge”. We have 
shown that the radial index is not tied exclusively to the 
transverse spatial profile of the beam but also has a phys¬ 
ical meaning. We hope that this opens up a new area of 
investigation and inspires some new potential technolog¬ 
ical prospects. 

We briefly outlined one such potential application, 
which is that the radial index may be naturally preserved 
under propagation in a graded-index fiber. This may 
have application to fiber-optic multiplexing in the classi¬ 
cal and quantum domains. 

In conclusion we conjecture that, despite the resources 
devoted to the study of the orbital angular momen¬ 
tum degree of freedom, the radial degree of freedom is 
vastly richer due to its mathematical complexity, and the 
plethora of questions that remain unanswered. 
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The function Z satisfies the wave equation 


dt 



Z (x, t) = 0. 


(42) 


Where is the full laplacian. The potential Z 
can, itself, be recast in a more tractable form: Z = 
(ei, 62 , e 3 )x(a;, t), where x is a field and the e’s may be 
chosen to suit the situation at hand (i.e. the desired sym¬ 
metries of the system to be studied), for beams propa¬ 


gating in a straight line Z = (0,0, l)x{x,t), is common. 

All solutions to the wave equation - Eq. (42) - may be 
expressed as a superposition of plane waves, so 
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X{x, t) 


J dkN{k) 


(43) 


X 








Where N{k) represents a normalization factor. Then the 
RS vector can be written as 


F {x, t) = J dk e{k) (44) 

Where e is a fc-dependent complex polarization vector 
which is determined by the choice of the vector part of 
the potential Z and includes the normalization factor. 
This represents a gauge freedom. Note that the polariza¬ 
tion vector factors completely from the rest of the expres¬ 
sion. It can thus be considered a part of the transforma¬ 
tion from momentum space to position space. This is a 
special feature of the RS formulation. The weight func¬ 
tions 'ijj~{k) and ip'^^k) are the positive and negative 
frequency components of the RS vector in momentum- 
space, the choice of which completely define the physical 
degrees of freedom of the EM field in momentum space. 

In our treatment we wish for a fully well-formed con¬ 
struction, so we will have as our objective to “end up” in 
momentum-space, since it has been long-known that the 
definition of the photon wave-function in position space 
is problematic. This is due to the fact that a position 
operator does not exist for photons (eigenstates of mo¬ 
mentum for photons - plane waves - though non-physical 
can be mathematically tractable, but the idea of a posi¬ 
tion eigenstate for photons is non-nonsensical even at the 
conceptual level as photons are inherently de-localized). 
A probabilistic position representation can generated for 
photons but we will not take this route. For more infor¬ 
mation on this and connected topics please see Refs. 00. 

However, for reasons of clarity we will start in posi¬ 
tion space and then switch to momentum. The time- 
dependent Schrodinger equation is 


ih —'ii {x,t) = H ^ (x, t). 


(45) 
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Written in index summation notation the tensor s = Sijk 
is both equivalent to the Levi-Civita symbol Cijk times 
a factor of —i, and the spin-1 matrices in quantum me¬ 
chanics (when these matrices are acting on the Carte¬ 
sian vector components of the wave-function and not the 
eigenstates of s^). Since it is the case that axb = Cijkajbk 
and a ■ b = Qibi, it is also the case that 


V X 4' = -if - V 4'. (46) 

Where p = —ihV is the very well-known momentum op¬ 
erator. Slightly less well-known is the “helicity” operator 
in optics and particle physics, which is the sign of the 
projection of the angular momentum on the momentum 


A = Sgn 

M ■ p 

5 


= Sgn 

(^—ihf X V -1- • 

i-iHV) 


1 . 

= —s ■ p. 
Pt 


(47) 

(48) 

(49) 


Where the momentum and angular momentum opera¬ 
tors are substituted for their definitions on the second 
line (the later is composed of an orbital and spin part), 

Pt = \jpx + Py + PzJ the simplification on the last 

line is due to the fact that (r x V) • V = 0. So we 
can see directly by comparing the above to Eq.([4^ that 
the helicity operator (up to constants) is the Hamilto¬ 
nian. The helicity operator generates the duality rota¬ 
tion and is associated with that symmetry in free-space. 
Other than the fact that it serves as the Hamiltonian 
for a Schrodinger equation with the RS vectors as wave- 
functions, the properties of A are not strictly relevant to 
our work here, nonetheless the interested reader is re¬ 
directed to Refs. [26] [27] . 

It will be convenient to find the eigenfunctions of the 
helicity (and thus the energy eigenstates) in a beam-like 
basis. We return to Eq.([4^ and utilize the following 
expansion in cylindrical coordinates 


Where FI is the Hamiltonian. One immediately recog¬ 
nizes the similarity between this and Eq. (39), the first 
Maxwell’s equation for the RS vector. We can thus 
equate the Hamiltonian to the appropriate part of that 
equation = —icV x 4^ with the RS vector becoming 
the photonic wave-function. 

In order to proceed from here we will rewrite the curl 
operator. Take a vector s, the components of which are 
the spin-1 matrices 




^ (50) 


Where k^ and kt are the polar and radial coordinates 
in momentum space, respectively. Note that the radial 
coordinate in momentum space is also the transverse mo¬ 
mentum. The functions Jm are Bessel functions of the 
first kind. This affects a change of basis from plane waves 
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to Bessel waves, which also forms a complete set. Now aFl 


ktkz 


The wave-functions in Eq.(53) form a basis of 


Eq.(43l becomes 


Xix,t) f dk x'L{k)'(p'^ik). (51) 

m,cr 

Where a can take values ±1 and the Xm’s are defined as 


X^{x,k) = (52) 

We then obtain via choosing the beam-like form of the 
potential {Z = (1,0,0)x) and by direct calculation for a 
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JmiktX ). 


solutions to the photonic Schrodinger’s equation. It is im¬ 
portant to note that, since the helicity takes the place of 
the Hamiltonian, there are solutions to the Schrddinger’s 
equation which at first glance might seem to have neg¬ 
ative energy. These solutions (those of left-circularly- 
polarized photons) do not represent anti-particles as pho¬ 
tons have no anti-particles. They should be interpreted 
merely as opposite helicity photons. These beam-like so¬ 
lutions are in fact the Bessel beams, which have come to 
be objects of interest due to their diffraction-free proper¬ 
ties [IH]. However, much like plane-waves, these beams 
are not physical as they turn out to have infinite en¬ 
ergy. Typically these beams are made physical without 
breaking the assumption of monochormicity by convolv¬ 
ing them with a Gaussian function (the Bessel-Gauss 
beams). However these beams are no longer exact so¬ 
lutions to the Maxwell’s equations. However by drop¬ 
ping the assumption of a monochromatic beam (also usu¬ 
ally implicit in the paraxial wave-equation) we can use 
the Bessel-beam basis defined in Eq.([5^ and continue to 
have exact solutions to the photon wave-function which 
are valid under all circumstances. 


The number a represents whether the photon is right or 
left circularly polarized. It is the case that = 












